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Abstract 

The Short GAmma Ray Front Air Cherenkov Experiment (SGARFACE) uses the Whipple 10 m telescope to search for bursts 
of 7 rays. SGARFACE is sensitive to bursts with duration from a few ns to ~20 /is and with 7-ray energy above 100 MeV. 
SGARFACE began operating in March 2003 and has collected 2.2 million events during an exposure time of 2267 hours. A 
search for bursts of 7 rays from explosions of primordial black holes (PBH) was carried out. A Hagedorn-type PBH explosion 
is predicted to be visible within 60 pc of Earth. Background events were caused by cosmic rays and by atmospheric phenomena 
and their rejection was accomplished to a large extent using the time-resolved images. No unambiguous detection of bursts of 7 
rays could be made as the remaining background events mimic the expected shape and time development of bursts. Upper limits 
on the PBH explosion rate were derived from the SGARFACE data and are compared to previous and future experiments. We 
note that a future array of large wide-field air-Cherenkov telescopes equipped with a SGARFACE-like trigger would be able to 
operate background-free with a 20 to 30 times higher sensitivity for PBH explosions. 
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1. Introduction 

SGARFACE is designed to detect short bursts of 
7 rays with energies above 100 MeV using the air- 
Cherenkov technique The design allows detection 
of bursts lasting up to ~20 /is and has a directional ac- 
curacy of up to a few arc-minutes [2J. Bursts of this 
type could be produced by the evaporation of primor- 
dial black holes (PBH), from emission associated with 
giant radio pulses from pulsars J3), or by a high-energy 
component accompanying very short gamma-ray bursts. 
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Primordial black holes would have been the first ob- 
jects formed by the gravitational collapse of inhomo- 
geneities in the early universe 0415161 . Inhomogeneities 
may have resulted from density fluctuations [7| present 
after inflation [8 1, the presence of cosmic strings 1191 101 , 
collapse of domain walls IfTTl . or from bubble collisions 
during a phase transition [ 12 1. Consequently, even with- 
out a detection, limits on the PBH density provide im- 
portant insight into many aspects of the physics of the 
early universe |[T3l , including: cosmological phase tran- 
sitions 1 14 1 and references therein, structure formation 
0151161 . mass clustering 0171 . and magnetic monopoles 

Da. 

Three effects have been used to search for PBHs: 
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(1) Gravitational effects: Black holes of mass > 10 -8 
Mq (2 x 10 25 g) may be detectable by their gravita- 
tional effects in microlensing observations 1191201 or 
from mass clustering measured with the Lya forest |[T5l . 
A significant number of heavier PBHs may have formed 
binary systems 1211 and may be detected from the gravi- 
tational radiation emitted during coalescence. An upper 
limit on the BH density in the mass range from 0.2 to 
1.0 Mq has been reported by the LIGO collaboration 

(2) Hawking radiation: Black holes of mass less than 
10 14 g may be detected by their emitted Hawking ra- 
diation [23 1. Black holes of mass, M, emit real parti- 
cles with a black-body energy spectrum at temperature 
T = 1.06(10 13 g/M) GeV E3). As mass accretion by 
PBHs is thought to be negligible [25|, PBHs will evapo- 
rate completely in 10 64 (Af/M Q ) 3 years, where Mq rs 
2 x 10 33 g is the Solar mass. At the present time, black 
holes with an initial mass of 4.5 x 10 14 g [26 1 would 
be reaching the final phase of their evaporation which 
might be explosive. PBHs that have already evaporated, 
would have left an imprint on the cosmic particle back- 
grounds and set stringent limits on the PBH density; a 
discussion is given in Sect. 11.11 Lastly, PBHs may also 
be detectable as transient gamma-ray sources in the so- 
lar system H27I28I . 

(3) Accretion effects: The accretion onto light black 
holes may produce distinct observable radiation in the 
relatively dense environments of interstellar space, in 
binary systems, or near planets [29 1 . 

The types of particles emitted by Hawking radiation 
and their energy spectra are well understood below the 
QCD confinement temperature of T c w 175+^ 8 MeV 
ll30l . Above this temperature, it is unknown whether 
individual elementary particles will be emitted in the 
strong gravitational field, or if a phase transition of 
the surrounding vacuum occurs and a thermodynamic 
model of the emission is valid. If the standard parti- 
cle physics model (SM) is valid at all temperatures, the 
emission spectra can be determined from the quark- 
gluon jet fragmentation functions convolved with the 
Hawking emission function [24 26 1. This would result 
in a slow evaporation at high temperatures: the final ex- 
plosion would last about 1 s with 7-rays of energies 
above 400 GeV being emitted 1311321 . 22% of the total 
energy would be emitted in 7 rays with a flux spectrum 
peaking at ~ 100 MeV and falling off as E~ 3 at higher 
energies l26l . 

On the other hand, in thermodynamic models of 
the QCD phase transition point, the large number of 
hadronic resonances implies a large phase space and 
therefore the final evaporation becomes extremely 



rapid, resulting in an explosion. The most extreme ther- 
modynamic model is that of [33 34|, referred to as the 
Hagedorn model (HM), where the density of states, p s , 
increases exponentially with mass: p s oc m~ 5 / 2 e m / Tc . 
In this picture, the ground state of the vacuum contains 
elementary and composite particles, which are evapo- 
rated by the gravitational field of the PBH. During the 
final explosion, 10-30% of the mass at T c is converted 
into gamma-rays between 100 MeV and 1 GeV; the 
explosion lasts ^100 ns [27]. 

Instead of completely evaporating, PBHs may leave 
behind Planck mass relics, see 1T31 and references 
therein. Because of their small mass, Planck relics do 
not affect 7-ray emission, but may possibly contribute 
to dark matter [35 36 1. If extra dimensions exist below 
the currently measured scale of gravity, ^10 /im (37], 
black hole evaporation would release less detectable 
energy Il38l39ll . 

For SGARFACE, we consider the detection of PBHs 
via their final-stage 7-ray radiation in the Hagedorn 
model. The nominal burst is assumed to come from a 
black hole of mass 6.5 x 10 13 g, corresponding to T = 
160 MeV. It is assumed that 20% of the mass-energy, 
or 7.5 x 10 45 eV, is emitted in the form of 1 GeV 7 
rays. Realistically, the 7-ray spectrum would not be a 
(^-function, but a product of the direct 7 emission with 
the decays from hadrons. It is straightforward to scale 
the upper limits on the PBH explosion rate presented 
in Section[3]with the total burst energy and the average 
energy of 7 rays, see Section l2~6l 

1.1. Summary of PBH Searches 

Limits on the PBH explosion rate have been placed 
by direct searches for the final stage emission and re- 
quire a particular particle physics model. Limits on 
the PBH density that have been derived from cosmic- 
particle background spectra require an assumed initial 
PBH mass spectrum. To compare the results, assump- 
tions need to be made on PBH clustering in the galac- 
tic halo [31 40 17] and, in the case of charged particle 
backgrounds, the enhancement due to confinement by 
the Galactic magnetic field 11281 . 

Air-Cherenkov telescopes (ACT), due to their large 
collection area, but small field of view of ^4°, are well 
suited to search for individual explosions. The sensitiv- 
ity to detect standard model PBH explosions is limited 
to distances of 0.4 pc l32l . and to about 200 pc for the 
HM bursts. Upper limits on the explosion rate are shown 
in Fig. Q] where Hagedorn burst models have been as- 
sumed at time scales shorter than 1 ps and SM bursts at 
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longer time scales. The first search for HM bursts was 
carried out by 041 142143 L who arrived at an upper limit 
for bursts of 150 ns duration. The results of their two 
experiments were carefully reanalyzed using air-shower 
simulations to characterize the trigger threshold across 
the field of view. The revised value of their more sen- 
sitive, short-baseline, experiment is shown in Fig. Q] 

The first search for SM bursts of duration between 
0.01 s and 1 s and energy above 7 TeV was carried 
out by [44 1 with a pair of scintillation particle detectors 
separated by 250 km. Several groups, using single or 
closely spaced detectors, have set limits at other energies 
and duration 0451461471481321491 . Some of the limits 
shown in Fig.Q]were standardized to the same SM burst 
by OT1 . Satellites, though they have a small collection 
area, are useful in searching for PBH explosions because 
they are essentially background-free and have a large 
field of view. EGRET data was searched for multiple 
7-ray tracks within its trigger window of 600±100 ns 
ll50l . Using BATSE data, a possible association of PBH 
explosions was made with some short (6-200 ms) X-ray 
bursts on the basis of the spectral evolution and spatial 
isotropy OP . 

The evaporation of PBHs over the lifetime of the uni- 
verse would have left an imprint on the observed cosmic 
particle spectra: 7, e ± , p, p, and t/P"lThe underlying as- 
sumption in placing limits on the PBH density using par- 
ticle backgrounds is that the initial PBH mass distribu- 
tion arose from scale-invariant density fluctuations l27l . 
Using the measured anisotropy in the EGRET 7-ray 
background above 100 MeV, [40 1 determined the rate 
of PBH explosions to be less than 0.07-0.42 pc~ 3 yr _1 , 
the confidence level of this upper limit is not given. Of 
the charged particles, p are the most effective in setting 
limits on the PBH density M28I52I . however additional 
uncertainties on the leakage time out of the Galaxy are 
introduced. The limit on the PBH explosion rate from 
the p data is (0.01 1-0.067) pc -3 yr _1 , where p are pro- 
duced by jet-fragmentation. Here, the systematic error 
due to uncertainties in the halo enhancement may be 
as large as a factor of two ll52l . The most recent and 
stringent upper limit comes from [53 1, who use a 3D 
Monte Carlo code to model the p propagation in the lo- 
cal Galactic environment and compare it to the BESS 
1995 p data (54]. We scale their result to the 99% CL as 
3.4 x 10~ 2 pc -3 yr _1 , assuming Poisson statistics. An- 
other method to detect the evaporation of PBH relies on 
the neutrino background. The flux emitted in neutrinos 
is at least 2 to 3 times larger than the 7-ray flux lf55l . 



1 We use the following constants: Hq=11 km/s/kpc, 13.6 X 10 9 yr 
as the age of the universe, and a critical density of 1 X 10~ 29 g/cm3. 



but so far only the low-energy ^-background has been 
used to constrain the PBH mass-distribution [56|. 

Another idea to detect PBH explosions was put for- 
ward by ||571 , who suggested that a HM burst produces 
a sphere of simultaneously emitted charged particles, 
but see the arguments by |58| and |59|. These generate 
radio and optical emission by interaction with the am- 
bient magnetic field. Searches for optical pulses have 
been carried out by [60 42 ] resulting in an upper limit 
of 0.3 pc~ 3 yr -1 at the 99% CL. A dedicated search 
for coincident optical flashes between two optical tele- 
scopes with duration between 1 and 30 ms by lETTl set 
an upper limit at 0.028 pc -3 yr -1 . Searches for radio 
pulses from PBH evaporation are summarized in BTI 
; the best limit being < 7 x 10~ 10 pc~ 3 yr -1 . How- 
ever, as the physical interaction model of [57 1 is highly 
speculative, these upper limits should not be regarded 
as definite. 
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Fig. 1. Limits on the density of PBH explosions. Direct-detection 
limits are at 99% CL, except for 1511 . and use fl m = 0.06. The 
limits from particle background measurements have possibly large 
systematic uncertainties associated with them. 



1.2. Description of SGARFACE 

SGARFACE J2] is now operating as part of the Whip- 
ple 10 m telescope, an imaging air-Cherenkov tele- 
scope (ACT) located in southern Arizona [62 63 1. The 
Whipple 10 m telescope is located at 31.6804° latitude, 
110.8790° W longitude, and 2312 m a.s.l. 

The camera of the Whipple 10 m telescope consists 
of 379 close-packed PMTs covering a 2.4° field of view. 
The PMT signals for SGARFACE are split off the 10 m 
signal cables via passive couplers. To reduce cost and 
complexity of the SGARFACE system, clusters of seven 
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nearest-neighbor pixels are summed into 55 channels, 
each viewing ^0.36° of the sky. This pixelation is suf- 
ficient to image the extended Cherenkov images pro- 
duced by bursts of 7 rays. The channels are sampled by 
flash analog-to-digital (FADC) converters at an inter- 
val of 20 ns and with memory depth of 35.04 /is (1752 
samples). A trigger, described below, causes the FADC 
buffers to be stored to disk. The FADC input was de- 
signed with a time-constant of about 25 ns to preserve 
the accumulated charge between samples. SGARFACE 
achieves background rejection of cosmic rays and at- 
mospheric phenomena by the time-resolved imaging of 
the Cherenkov wave fronts. 

The air-Cherenkov technique has a detection area of 
~ 1 x 10 5 m 2 for single TeV 7-ray primaries l64l . 
At sub-GeV energies, the Cherenkov emission from 
one 7 ray is not detectable, but a wavefront of near- 
simultaneous sub-GeV 7-rays produces a detectable sig- 
nal (Tj. The time dispersion of Cherenkov photons from 
an instantaneous burst is 10 ns (35 ns) on-axis (3° off- 
axis) due to the shower geometry, while a 6.5 ns time- 
spread is introduced by the geometry of the 10 m tele- 
scope. The angular distribution of Cherenkov photons 
of these low-energy 7-ray wavefronts is smooth and 
broad, as shown in Fig.[2]for a burst of 1 GeV 7 rays at 
zenith. The widths of the angular distributions, shown 
in the lower part of Fig. are a combination of Comp- 
ton and Coulomb scattering, and the geomagnetic bend- 
ing of e . Taking into account the number of radia- 
tion lengths available for electrons until they reach the 
threshold for Cherenkov emission in the atmosphere, 
the total multiple Coulomb scattering angle is roughly 
proportional to p^ 1 ^ 2 , where p is the initial electron 
momentum [65|. Similarly, the width of the photon an- 
gular distribution arising from Compton scattering is 
proportional to p~ 1//2 , where p is the photon momen- 
tum ll66l pg. 354-358]. As the primary 7-ray energy de- 
creases towards the Cherenkov threshold, relatively less 
Cherenkov radiation is produced from the total avail- 
able energy; the threshold for e^ 1 is 80 MeV at 20 km 
altitude, the height of shower maximum for 100 MeV 
7 rays. Therefore, the Cherenkov-light density on the 
ground depends strongly on the primary energy. For a 
fixed total energy, Figure [3] shows the Cherenkov pho- 
ton yield with primary energy relative to 1 GeV 7 rays. 

The SGARFACE trigger consists of two systems J2j. 
A multi-time scale discriminator (MTSD) and a pattern- 
sensitive unit. The MTSD integrates the signal and trig- 
gers on six time scales: 60 ns, 180 ns, 540 ns, 1620 ns, 
4,860 ns, and 14,580 ns. The trigger threshold require- 
ment is tested at three equally spaced time intervals dur- 
ing the trigger window, requiring pulses to be at least 
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Fig. 2. Angular distribution of Cherenkov light received at 2300 m 
altitude from a plane wave front of 1 GeV 7 rays at zenith with 
fluence of 4 7/m 2 . The x-axis is perpendicular to the Earth's mag- 
netic field. The solid line shows the photon distributions along the 
x-axis; the dashed line along the y-axis. 
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Fig. 3. Variation of the Cherenkov photon yield divided by the pri- 
mary energy, Eq, versus Eq, corresponding to the yield of Cherenkov 
photons for a burst with fixed total energy. The yield has been nor- 
malized relative to Eo = 1 GeV. The yield decreases as the primary 
7-ray energy approaches the threshold for Cherenkov emission. 
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as long as the trigger window. The thresholds were set 
conservatively above the night-sky background and oc- 
casionally adjusted to reflect changes in the telescope 
throughput. 

The pattern-sensitive coincidence unit determines if 
a sufficient number of nearest neighbor channels have 
triggered the MTSD. The majority of observation time 
was taken with a 7-fold nearest-neighbor coincidence 
requirement. Once an event triggers, the 55 FADC 
stacks are read out after a 16.28 fis delay. This data is 
stored together with information on telescope pointing, 
PMT voltage settings, and trigger information. 

2. Data Analysis 

2.1. Overview 

Since March 2003, SGARFACE has recorded about 
2.2 million events during routine observations of the 
Whipple 10 m telescope. Observations with the 10 m 
telescope are usually only carried out when the weather 
conditions are judged acceptable for TeV astronomy ob- 
servations. The mean data rate was 0.18 Hz. The dead- 
time of the data-acquisition (DAQ) system is 0.316 sec 
per event, and comprised 6% of the total operational 
time. The lifetime accumulation is shown in Fig. [4] 
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Fig. 4. Lifetime accumulation occurs when the following conditions 
were met: trigger rate greater than 3 events/min, telescope pointing 
above 20° in elevation, telescope pointing varied by less than 4° 
root-mean square during the one minute interval. The last condition 
is not met when the telescope is being slewed to a new source. 

The exposure of the sky during the six years of op- 
eration is shown in Figures [5] in equatorial, galactic, 
and horizon coordinates, respectively. Long exposures 
at prominent TeV gamma-ray sources are at the Crab 
Nebula, Mrk 421, Mrk 501, the Galactic Center, and 
H 1426+428. Near each source is an exposure off- 
set in Right Ascension (RA) by ±30 min that arises 
from Whipple OFF-source observations. The horizon- 



tal bands in the equatorial exposure map in Fig. [5] stem 
from the nightly 10 minute zenith runs when the tele- 
scope is pointed at or near the zenith. 
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Fig. 5. Exposure map (minutes) in equatorial (top), galac- 
tic (middle), and horizon (bottom) coordinates. For visual clar- 
ity, bins are suppressed with less than 5 min exposure. 
Some prominent points of observation are the Crab Nebula 
(RA 05h34m31.9s DEC+22°00m52.1s, J2000), Mrk 421 (RA 
Ilh04m27.31s DEC+38° 12m31.8s), and the Galactic Center (RA 
17h45m40.0s, DEC-29°00m28.1s). The Galactic Center is near 360° 
longitude 



2.2. Event Analysis 

For each event, the topology and time development 
of the images is parameterized on five time scales. As 
in the hardware algorithm, the analysis sums the FADC 
traces with a top-hat filter function of width 1, 3, 9, 27, 
and 81 samples. The analysis then proceeds on each 
time scale as follows: 

Fluctuations about the pedestal value are produced 
by night-sky background fluctuations, electronic noise 
from telescope-drive motors, and intrinsic electronics 
noise of the DAQ. In the analysis, the mean pedestal 



5 



value is subtracted from the trace. The standard devia- 
tion of the distribution of the fluctuations, referred to as 
RMS in the following, is used to set the threshold that 
channels must cross before a particular sample is ac- 
cepted as part of the image. Details of the pedestal sta- 
bility are given in Sect l6.2l The triggering sample in the 
FADC stack is located at least 814 samples before the 
end of the trace. Therefore, the pedestal of each channel 
is determined from the first 750 samples with a value 
less than 40 dc multiplied by the width of the top-hat 
filter. The longest trigger time scale of three times 243 
samples, or 14.58 /is, is not used in the analysis because 
the pedestal value cannot be estimated accurately from 
the initial signal-free portion of the FADC trace. The 
loss in fiuence sensitivity by excluding this time scale 
is minimal. 

Next, FADC samples that contain little or no signal 
are suppressed followed by a moment-analysis to char- 
acterize the image shapes. The cleaning procedure pro- 
ceeds in two steps: 1) Any re-binned samples with sig- 
nal in excess of three times the pedestal RMS value are 
retained. 2) Samples adjacent in time and space are re- 
tained if their signal is greater than twice their pedestal 
RMS value. All other samples are ignored in the anal- 
ysis. 

The analysis evaluates the trigger requirement on five 
time scales with a fixed set of thresholds for the entire 
6-year data set. The thresholds were chosen just above 
the highest hardware levels used during the majority of 
observations: the thresholds on the trigger time scales 
are 18 dc (60 ns), 53 dc (180 ns), 119 dc (540 ns), 325 
dc (1620 ns), and 874 dc (4860 ns). 

Details of the image and timing parameterization are 
presented next. The parameterization proceeds along the 
following steps: 

1) The duration of the pulse is measured by the 1/e- 
width of the summed FADC channels. 

2) The event is integrated for the duration of the pulse. 
The integrated image is characterized by the centroid 
of the light distribution, the width and length, and the 
rotation angle, a, of the major axis away from the cam- 
era center, see Fig|6] These parameters are calculated 
from the moments of the light distribution |67 68). In 
addition, the peak values of the seven brightest pixels, 
max\..T, are recorded. The contrast ratio is measured 
from the ratio of the mean brightness of pixels located at 
distances less than 0.4° from the centroid, versus those 
at distances greater than 1°. 

3) The speed of the shower core along the image ma- 
jor and minor axes is measured from the instantaneous 
centroids of the ten brightest pixels. 




Fig. 6. Illustration of Cherenkov radiation imaged by SGARFACE. 
Pixels are shown as thin outlines and the amount of light is rep- 
resented by filled circles. The shower image is parameterized by 
a two-dimensional Gaussian distribution with a measured width, 
length, centroid, and orientation angle, a, as indicated. 

2.3. Monte-Carlo Simulations 

The sensitivity of SGARFACE is determined from 
Monte-Carlo simulations. We use the KASCADE J69) 
air-shower simulation program, version 7.3, together 
with GrlSU(tah) [70| to simulate the air-Cherenkov pro- 
duction and telescope response. The calibration of the 
light-throughput of the telescope and the gain of the 
PMTs is described in Section |6*TI 

The simulation of the detector electronics includes 
the single photo-electron pulse shape received by 
SGARFACE and is shown in Fig. [7] The long tail in the 
pulse shape is caused by the implementation of the AC 
coupling of the signal cables and their back-termination 
at the camera. Due to this signal shape, more charge 
is accumulated on longer integration time scales. This 
reduces the effective trigger threshold and produces an 
increase in sensitivity nearly proportional to V time. 
This was unforeseen in the design of the trigger and 
helps the performance of the MTSD. The trigger thresh- 
old at the shortest time scales corresponds to about 12 
Cherenkov photons per m 2 per pixel per 60 ns. 

2.4. Background Rejection 

Data from 2003/4 was used to understand the back- 
ground events and, in conjunction with simulated bursts 
of 7 rays, to develop a background-rejection method. 
Several sources of background are seen in the data and 
can be separated into atmospheric (cosmic rays, light- 
ning) and man-made phenomena. 

The largest number of background events are cosmic 
rays where the Cherenkov radiation is seen at various 
off-axis directions resulting in event durations of up to 
100 ns. SGARFACE also detects events lasting up to 
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Fig. 7. The measured pulse shape of single photo electrons at the 
input of the FADCs. 



^450 ns where a central bright spot moves through the 
field of view. These events are compatible with fluores- 
cence emission from cosmic rays. 

From comments by the 10 m operators, we deter- 
mined that lightning produces another type of back- 
ground seen occasionally in the data. These events have 
a patchy structure and are longer than 16 /us, compati- 
ble with the measurements by [71] that show a duration 
between 25 [is and 2 ms. Another type of background 
that we considered is light from meteor trails. Though 
the light in the trails has a duration long enough to trig- 
ger the MTS, ^0.3 s [72], the lateral angular extend of 
< <0. 1 1731 is not large enough to trigger the pattern- 
sensitive unit with a 7-fold coincidence requirement. 

The 10 m observing log contains the rating of weather 
conditions by the observers. We did not use these rat- 
ings in cutting run-time a-priori for two reasons: One, 
the ratings are subjective and two, the GPS time in- 
formation was not recorded properly in the 10 m run- 
header files beginning in 2007. The later would have 
made a cross-correlation of the weather information 
with SGARFACE events a very labor intensive project. 
Instead, we decided to use all available data and only in- 
vestigate the observing conditions for events that passed 
our selection criteria. The majority of 10 m observa- 
tions are taken during good or very good night-sky con- 
ditions, though occasionally the telescope was operated 
for short periods of time when the weather was bad. By 
including all available data, the integrated lifetime of 
SGARFACE may be too large large by as much as 10%. 
We do not apply a compensating lifetime reduction in 
the analysis. 

Background that was identified as man-made includes 
the light produced by airport/airplane beacons and is 
detected either directly or from the reflection by haze 
or clouds in the night sky. Directly visible from the 
Whipple 10 m telescope is the beacon light at the No- 
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Elevation 


>20° 


Trigger 


7-fold 


maxi 


<230 dc 


Pulse stop 


<1750 - W 



Table 1 

Preliminary cuts applied to the entire data set. Here, W refers to 
the width (in 20 ns samples) of the top-hat filter function. 

gales airport. We have identified some events as arising 
from these beacons based on the repetition rate. A third 
source of background is RF noise picked up by the data 
acquisition system. This background typically affects 
only a small fraction of the signal channels, but can 
cause triggers on long time scales. To eliminate these 
known sources of background, events are required to 
pass the "preliminary" cuts shown in Tab. Q] 

Following the application of these preliminary cuts 
to the entire data set, 100,221 events remained. A more 
selective set of cuts was then developed to reject as 
much background as possible while retaining all simu- 
lated bursts. Bursts were simulated at discrete energies 
between 100 MeV and 10 GeV, off-axis arrival direction 
of up to 4°, and fluence from 0. 175 to 32 7-rays/m 2 . The 
simulated telescope pointing directions are as shown in 
Fig. El 

For any particular telescope pointing and burst prop- 
erty, the measured burst parameters occupy a narrow 
space compared to the data. Fig.[8]compares the param- 
eter distributions of data taken in 2003 with simulated 
bursts of 7 rays for one set of simulation inputs de- 
scribed in the caption. The parameter distributions do 
not change significantly with the width of the top-hat 
filter or with the burst duration. In principle, this separa- 
tion in parameter space between bursts and background 
could be used to achieve a very efficient background 
rejection. 

The fluence and location of the burst in the field of 
view and, at lower elevations the azimuth dependence, 
create a five-dimensional parameter space for which 
specific cuts would need to be developed. We do not 
consider it practical to do so at this time. An alterna- 
tive, which we also do not implement because of its 
dramatic effect on the phase space to detect PBH explo- 
sions, would be to restrict the data set to a very limited 
field of view and telescope pointing. Instead, we devel- 
oped a very simple set of cuts, given in Table that 
does not reject any simulated bursts of gamma rays. We 
note that with these simple cuts we expect some back- 
ground events to pass the cuts and these events will be 
investigated a-posteriori. The cut on width with tele- 
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Parameter 


Cut 


Length 


<0.65° 


Distance 


0.47°... 0.58° 


Contrast 


<0.67 


Centroid speed 




time scale 1 


<0.0276°/ns 


time scale 3 


<0.0109°/ns 



Table 2 

Cuts on the parameters that retain all simulated bursts. 

scope elevation is shown in Fig. [9] and compared to the 
width seen for data. This set of cuts accepts all simu- 
lated bursts of 7 rays and retains 36,457 events out of 
the entire data set. 

The variation of the mean centroid position, distance, 
with off-axis arrival direction, is shown in Fig.[l0] The 
burst direction can be determined from the distance pa- 
rameter, if the distance is less than 0.25°. At larger 
pointing offsets, the field of view limits the positional 
reconstruction ability. 

2.5. Detector Sensitivity 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Width Length (°) 




\0~* 10~ z 

Centroid speed (°/ns) 



0.2 0.4 0.6 0.8 1 

Contrast ratio 



The fiuence sensitivity is determined as a function 
of burst duration, r, primary 7-ray energy, Eo, the tele- 
scope elevation, El, azimuth, AZ, and the burst posi- 
tion within the field of view, 9. We make the assump- 
tion that the sensitivity can be expressed in terms of 
four independent functions 



Fmin [El, AZ, T,Eq,0] — 



P 



g[E ] x h[El,Az] x k[9] ' 
(1) 

where F m i n [t] is determined after preliminary cuts and 
at nominal conditions: telescope pointing at zenith, So. 
= 1 GeV, and the burst is centered in the field of view, 
6=0. The upper limit on the detectable fiuence, F max , 
is due to the limited dynamic range of the FADCs, and 
it is determined from simulations. The lower and upper 
limits on the fiuence are determined by requiring a 90% 
trigger probability. The contribution of each trigger time 
scale to the fiuence sensitivity is shown in Fig.[TTl where 
only trigger time scales are shown that enhance the sen- 
sitivity at longer burst durations. Already at the shortest, 
60 ns, time scale, the fiuence threshold increases nearly 
at the noise limit, proportional to V time, as can be seen 
by comparison with the thin dotted line in Fig.Qj] The 
systematic uncertainty in the fiuence is due to the en- 
ergy calibration of the detector and is estimated to be 
15%. 



Fig. 8. Normalized parameter distributions of data taken during 2003 
(dashed) after application of preliminary cuts and at elevations above 
80° . The hatched areas show the parameter range of simulated bursts 
of 1 GeV 7 rays arriving on-axis at zenith with a burst fiuence of 2 
7/m 2 . Image parameters were calculated using a one-sample wide 
top-hat filter. The large excess of events in the histogram of contrast 
ratios ~0.8 is due to flashes of light from a Nitrogen pulser used 
in the nightly calibration procedure of the 10 m telescope. 



% 0.6- 

3 



W 0.5- 

E : 

0.4F 
0.3p 
0.2P 
0.lh 

c 



20 



30 



40 



50 



60 70 80 90 
Observing elevation (°) 



Fig. 9. Measured image width vs. telescope elevation for events seen 
in the 2003 data. The shaded region shows the rejected parameter 
space. The excess of events in the upper part of the figure (width 
of ~0.58° 78° observing elevation) are not of astrophysical origin 
as they occurred while engineering work was performed on one of 
high-voltage power supplies. 
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°0 0.5 1 1.5 2 2.5 3 3.5 

Simulated burst direction (°) 

Fig. 10. Variation of the measured distance of the image centroid 
from the center of the field relative to the true burst arrival direction. 
The error bar's show the RMS spread of the distribution. 

To determine the functions g[Eo], h[El,Az], and 
k[G], we note that the fluence threshold is proportional 
to the value of the seventh brightest pixel, max7. We use 
the relative change of maxy with respect to Eq=1 GeV, 
El = 90°, Az=0°, and 6=0, to determine these functions. 
The variation with primary energy was already shown 
in Fig. |3j here it is formally identified with g[Eo]/Eo- 

Cherenkov images from bursts of 7 rays at 1 GeV are 
of 2 to 3° angular extend, as shown in Fig. [2] SGAR- 
FACE has a geometric field of view of 2.4° and is sen- 
sitive to bursts originating well outside of this range. As 
burst images are slightly oblong due to the Earth's mag- 
netic field, the sensitivity decreases faster for an offset 
parallel the direction of the Earth's magnetic field than 
perpendicular to it, see [2|. We ignore this directional 
dependence and use the average value to characterize 
the off-axis sensitivity; i.e. k[9] = k[9], only. The relative 
sensitivity with respect to the burst direction in the field 
of view is shown in Fig. Q~3] The sensitivity variation 
with telescope pointing, h[El, AZ], shown in Fig. [12] is 
due to shower geometry, Cherenkov-light absorption in 
the atmosphere, and the effect of the geomagnetic field. 

The functions, g, h, and k are also used in scaling 
the upper limit, F max . This is not strictly correct as 
the maximum fluence is limited by the brightest pixel, 
maxi, and not max7. However, the introduced error is 
estimated to be small, 5-10%, and the impact is negli- 
gible on the total volume for detection of PBH explo- 
sions. 

2.6. Detection Volume for PBH Explosions 

The volume to detect PBH explosions is dV" = 
(i?^ ax - i^Jdfi/3, where df2 = d<£sin[f?]d0 is 
the solid angle, and R m in,max are the detection lim- 
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Duration (ns) 

Fig. 11. Accepted fluence versus burst duration at various trigger 
time scales for a 1 GeV burst of 7 rays observed at zenith and 
centered in the field of view. Only the three trigger time-scales 
that increase the acceptance are shown. The solid lines outline the 
trigger regions corresponding to each timescale. The contribution of 
the 1.62 fis trigger timescale is shown by the small, dark shaded 
area at the far right. The dotted line shows the fluence that could 
be achieved if the trigger were at the noise level, i.e. increasing 
proportional to Vtime. 
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Fig. 12. Trigger sensitivity relative to zenith as a function of 
telescope azimuth and elevation. North is at 0° , East is at 90° AZ. 
The sensitivity to detect bursts of gamma rays is greatest along 
the direction of the Earth's magnetic field. At the location of the 
Whipple 10 m telescope, the Earth's magnetic field points at 170.2° 
azimuth and 59° elevation. Data was provided by the Geomagnetism 
Department at the National Geophysical Data Center. 

its on the distance compatible with F m i n ^ max . The 
distance limits to detect bursts at nominal conditions 
are shown in Fig. [14] The fluence, F, of a burst 
at distance, R, and initial particle number, N is 
F = N Q / (AttR 2 ), hence the volume scales with F as 
V = df2/(3(47r) 3 / 2 ) x {Nq/F) 3 / 2 . The initial number 
of particles can be written in terms of the total energy 
released in 7 rays, Et- No = Et/Eq. The detection 
volume for PBH explosions can then be written as 

, v df> E T xg[E }xh[El,Az}xk[9} 3/2 
3(4tt) 3 / 2 [ E ' 

x (F™„M~ 3/2 - F max [T]-V 2 ). (2) 

The total volume to detect bursts, integrated over a 6° 
field of view, is shown in Fig. [15] as a function of burst 
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Simulated burst direction (°) 

Fig. 13. Variation of the sensitivity to detect bursts arriving off-axis. 
The arrival direction is along a line inclined by 45° with respect to 
the x-axis. In the analysis, we limit the half-field of view to <3°. 

duration. The size of the field of view was chosen some- 
what conservatively so that the sensitivity at the outer 
edge is at least 20% compared to the center. The sys- 
tematic error of the volume is 23%, stemming from the 
uncertainty in the energy calibration of the detector. 





50 


o 


30 


mits 


20 


u 


15 


c 






10 


5 





0.05 0.1 0.5 1 5 10 

Duration [fis] 

Fig. 14. Maximum, Umax, (solid) and minimum, R m i n , (dashed) 
detectable distance for PBH explosions observed at the center of 
the field of view and telescope pointing at zenith. The total burst 
energy is 7.5 X 10 45 eV released in the form of 1 GeV 7 rays. 



3. Search for PBH Explosions 

Most events trigger on several time scales, as shown 
by the overlap of trigger regions in Fig. QT| In search- 
ing for PBH explosions, events are identified using only 
the parameters of the shortest triggered time scale. An 
exception is made for very short events of duration less 
than 60 ns, which only trigger on longer time scales 
due to the tail of the single pe-distribution. These short 
events are removed a posteriori by a special cut consist- 
ing of otax7> 100 dc and duration <60 ns. 

The duration of events passing all cuts is shown in 
Fig- HH Due to the simple parameters and loose cuts, a 



0.05 0.1 0.5 1 5 10 

Duration [ps] 

Fig. 15. Volume for detection of PBH explosions with telescope 
pointing at zenith. The total burst energy is 7.5 X 10 45 eV released 
in the form of 1 GeV 7 rays. 

number of background events pass the cuts. Instead of 
further refining the cuts at this point, we describe the 
most likely origin of the remaining events. The hypoth- 
esis that the events could only have originated from a 
burst of gamma rays is tested by investigating the im- 
age and time-development in more detail and consider- 
ing the circumstances under which the event occurred: 
weather, tracking/slewing, comments by observer in the 
nightly log. 

Three distinct populations are visible in Fig. [19] The 
largest number of events are very short and bright cos- 
mic rays of at least 5 TeV and >3° arrival direction 
with respect to the telescope pointing. A second detec- 
tor separated by at least a few hundred meter would be 
needed distinguish cosmic rays from bursts of 7 rays. 
A second population peaks near 300 ns, while a third 
group of events is seen to have durations longer 7 /xs. 

Next, details are given about the second population 
of events, those of duration between 120 ns and 500 ns 
in Fig. [19] The time-resolved analysis of the 87 events 
contained in this group reveals that a large fraction ex- 
hibit a characteristic spot of ^0.3° diameter moving 
through the field of view at ~ 5 x 10~ 3 °/ns. The three 
events that were detected at less than 45° elevation de- 
velop in a downwards direction in the field of view. 
Events detected near zenith appear to have random di- 
rections. If the detected radiation of these events is flu- 
orescence emission from a cosmic ray traveling perpen- 
dicular to the line of sight, then the total light collected 
by SGARFACE corresponds to primary energies be- 
tween - 2.5 - 25 x 10 17 eV d pg. 288-290]. Because 
the timing and image parameters were calculated using 
a moment analysis, it is difficult to distinguish the mo- 
tion of a compact spot from the image of a burst of 7 
rays seen largely off-axis. A more sophisticated analysis 
coupled with cuts on the timing gradient as a function 
of off-axis distance and image brightness would likely 
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have eliminated these events. After identifying by eye 
51 events that clearly show the characteristic spot and 
removing those from the data set, the event durations of 
the total data set, corresponding to 36,118 events, are 
shown in Fig. |20] Of the remaining events in the sec- 
ond population, 22 show a light distribution that peaks 
at the camera edge and slowly moves in time, making a 
burst origin unlikely. The other 14 events satisfy the cri- 
teria for burst images; that is they exhibit a broad image 
shape and uniform time development. Further investi- 
gation shows that 12 of these events were seen mainly 
off-axis, resulting in ambiguity as to the determination 
of their origin. It is possible that all events have were 
caused by cosmic -ray fluorescence emission with the 
detected variations being due to differences in energy 
and the viewing angle. The two remaining events oc- 
curred at 10:47 UTC, 20 April 2007 and at 08:03 UTC, 
23 April 2007. The events show an extremely broad 
angular light distribution, > 2°, and uniform time de- 
velopment with rise time of 20 ns, fall time of 40 ns, 
and total duration of ~140 ns. The time-integrated im- 
ages are shown in Fig. [16] The first event occurred dur- 
ing routine observations of Mrk 501 at 80.5° El and 
332° AZ, and the second event was recorded just be- 
fore the nightly calibration (zenith) run was carried out 
with the telescope pointed at 80° El and 0° AZ; cor- 
responding to a random direction in the sky. They oc- 
curred during very good night sky conditions and no de- 
tector malfunction could be identified. SWIFT detected 
a gamma-ray burst at 6:18 UTC, 20 April, but no burst 
was detected on 23 April. We note that atmospheric flu- 
orescence from a burst of X-rays could produce short 
events with a very broad angular distribution [75 1. How- 
ever, we do not claim to have detected bursts of 7 rays 
or of X-rays, because we cannot completely rule out 
atmospheric phenomena. 





-d.5 n.s 

X[°] 



Fig. 16. Two unusual events that exhibit a very broad light dis- 
tribution and appear to be centered in the field of view. The time 
development of the events was uniform and they lasted approxi- 
mately 140 ns. The fitted ellipsoid and the direction of the major 
axis are superimposed on the camera pixels. 



Fig. Q2] Twelve of these events occurred within 3 min- 
utes of each other on 10 November, 2007, with the tele- 
scope pointing South at 28° elevation during very bad 
night-sky conditions: distant lightning was noted in the 
observer log. Unfortunately, it was not possible to filter 
out such data automatically because, as mentioned in 
Sect. 12.41 the GPS time stamps in the 10 m database that 
are necessary to correlate the observer's weather-rating 
with individual SGARFACE events were missing since 
2007. 

During the previous night and under good night-sky 
conditions, six of the 24 events occurred within one 
minute of each other, while the telescope was slewing to 
a new observing position. These events were probably 
caused by RF noise generated by the telescope tracking 
system. Another two events occurred during engineer- 
ing operation while the telescope was not scheduled for 
astronomical observations. One event occurred during a 
run flagged as bad by the observers due to a high trigger 
rate in part of the camera. Another event, Fig. [18] shows 
patchy structure and is incompatible with the smooth 
distribution of the integrated Cherenkov light expected 
from a burst of 7 rays, shown in Fig. [2] In any case, 
these events are excluded because the 10 m was not per- 
forming astronomical observations. The remaining two 
events show a rising time-development past the initial 
drop-off from the peak and are longer than the 16 fxs 
FADC depth. These events passed the analysis cuts be- 
cause only the initial peak is parameterized while the 
subsequent time development is ignored. A possible ex- 
planation for these events is distant, high-altitude light- 
ning. In summary, no events detected at the longest du- 
rations can be explained solely by a burst of 7 rays. 

In conclusion, based upon our study of events iden- 
tified as cosmic -ray background, atmospheric phenom- 
ena, or events triggered by RF noise, we cannot exclude 
that no bursts of 7 rays have been detected and that the 
entire excess is due to background. Upper limits on the 
true PBH explosion rate, shown in Fig. [21] were cal- 
culated using Poisson statistics of the remaining back- 
ground events shown in Fig. [20] where an exposure time 
of 1502 hours was used. This exposure time arises from 
correcting the true lifetime of 2267 hours for sensitiv- 
ity changes with elevation and optical throughput effi- 
ciency. The upper limits can be scaled to different total 
and primary energies using Eq. [2] and Fig. [3] 



4. Discussion and Summary 



Durations longer than 7 /is were measured for 24 
events. The arrival direction of the events is shown in 



The direct search for PBH explosions with SGAR- 
FACE covered a broad range of 7-ray energies and burst 
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Fig. 17. Map showing the origin of the 24 events detected with 
durations longer than 7 fxs. 
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Fig. 18. An event of 9.5 fis duration that passes all cuts but exhibits 
patchy structure that is not expected from a burst of 7 rays. 
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Fig. 19. Histogram of the duration of the 36,170 events passing all 
cuts described in the text. 

durations compatible with a Hagedorn-type black hole 
evaporation. The experiment is collecting data since 
March 2003 and makes use of time-resolved images of 
Cherenkov radiation to distinguish the candidate bursts 
of 7 rays from the much more numerous cosmic -ray 
initiated showers. In the analysis, the detector was mod- 
eled in detail and Monte-Carlo simulations of bursts of 
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Fig. 20. Histogram of event durations after removal of events 
visually identified as having a compact spot moving through the 
field of view. It is likely that these removed events were produced 
by fluorescence emission of cosmic rays. 
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Fig. 21. 99% upper limits placed on the density of PBH explosions 
by SGARFACE and by other experiments. All direct-detection limits 
have been scaled to a total energy of 7.5 X 10 45 eV released 
in 1 GeV 7 rays. The background-free SGARFACE limits are 
shown with a thick dashed line. Indirectly determined upper limits 
from background measurements are not necessarily at the 99% 
CL level due to large systematic uncertainties in the galactic halo 
enhancement. 

7 rays were used to determine the astrophysical volume 
over which PBH explosions could be seen. Using only 
a single imaging telescope, this is the first time such a 
detailed search has been carried out and a high back- 
ground suppression of 98.4% was achieved. The ma- 
jority of the remaining cosmic -ray background events 
are indistinguishable from burst events with a single 
telescope. Figure [21] shows the upper limits that could 
have been achieved in the absence of the cosmic-ray 
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background, i.e. if a second, identical detector had been 
placed at a large distance and operated in coincidence. 
Originally, SGARFACE was designed when the cam- 
era on the 10 m camera had a 4.8° field of view. With 
this larger field of view an almost three-fold increase in 
PBH detection volume would have been achieved. 

The GLAST satellite is sensitive to bursts of 7 rays 
above 0.03 GeV energy. While the detector has a dead- 
time of 26 /.ts, multiple photons can be detected within 
the ~ 1 fxs coincidence gate width, enabling a direct 
search for bursts of 7 rays. In determining the fluence 
sensitivity of GLAST to sub-/is burst of 7-rays, we re- 
quire a minimum of 10 photons to be collected. Assum- 
ing that no events are detected, we show the 99% CL 
upper limits on the PBH explosion rate for a 2 year ex- 
posure in Fig. El Shown in Fig. |22] is the 99% upper 
limit over the entire energy range of GLAST, where the 
total burst energy was taken as before as 7.5 x 10 45 eV. 
GLAST is most sensitive to burst of 7 rays at ~100 
MeV due to the relation between sampled volume and 
7-ray energy at a fixed total burst energy. 

In comparing the upper limits shown in Fig. |2T| it 
should be kept in mind that different regions of space 
are sampled by the different detection techniques. For 
example, GLAST is limited in distance to 14 pc to detect 
our nominal burst of 7 rays, while SGARFACE probes 
PBH explosions out to 60 pc. Limits from cosmic par- 
ticle backgrounds measure the PBH distribution on the 
scale of the Galaxy, while relying on the assumed mass 
distribution function, galactic halo enhancement, and in 
some cases the cosmic -ray propagation model, to set 
limits on the current rate of PBH explosions. These dif- 
ferent techniques complement each other in searching 
for all possible signatures of PBH evaporation. 

In this paper we considered bursts of 7 rays to be pro- 
duced during the final, explosive evaporation of PBHs. 
In itself, a burst detected with SGARFACE would not 
be sufficient to claim the detection of a PBH explosion, 
as bursts might be produced by the emission of 7 rays 
contemporaneous with giant radio pulses 1 3 1 or perhaps 
via some other, as yet undiscovered, mechanism. The 
smoking gun to claim the detection of a Hagedorn-type 
PBH explosion would include measurements of the neu- 
trino and graviton emission, as well as radio and opti- 
cal observations B57I581 . but see l59l for a refutal on 
such emission. Future searches for bursts should strive 
towards such a multi-messenger approach. 




E„ [GeV] 

Fig. 22. Upper limits on the PBH explosion rate that GLAST can 
set in 2 years of observation versus energy of the emitted 7 rays, 
Eo. The total burst energy released in 7 rays is 7.5 X 10 45 eV. 

5. Outlook 

SGARFACE is background limited and an unam- 
biguous detection can only be made if the cosmic-ray 
background can be suppressed. The coincidences from 
a second telescope separated by at least a few hundred 
meters can be used to veto events caused by cosmic-ray 
initiated showers and other local phenomena. The VER- 
ITAS array of air-Cherenkov telescopes (76) is located 
7 km from the Whipple 10 m telescope. While SGAR- 
FACE has a lower fluence sensitivity than VERITAS 
to detect /is-bursts, the data from simultaneous obser- 
vations with the two instruments can be used to reject 
short cosmic-ray events. The analysis of these simulta- 
neous observations is forthcoming and will focus on a 
search for bursts of 7 rays from point sources, such as 
the Crab pulsar. 

The sensitivity of the air-Cherenkov technique to de- 
tect PBH explosions can be improved by implement- 
ing the imaging, multi-timescale trigger pioneered with 
SGARFACE, on an array of telescopes. With a tele- 
scope separation of ~1 km, a global array trigger system 
would operate background free and reduce the single- 
telescope trigger threshold. A larger mirror area would 
further increase the sensitivity, while a wide field of 
view increases the volume searched for PBH explosions. 
For example, if the array consists of telescopes with 
20 m diameter, 6° field of view, and improved noise 
characteristics, the volume over which PBH explosions 
can be detected increases by a factor of 20 to 30 over 
SGARFACE. 
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6. Appendix 

The accuracy to which the energy response of SGAR- 
FACE can be calibrated depends on the performance and 
stability of the optical throughput of the telescope, the 
photo-electron to digital-counts conversion ratio, and 
the uniformity of the camera response across the field 
of view. 

6.1. Calibration of the Light-Throughput 

The telescope response to Cherenkov light was cal- 
ibrated by direct measurement of all components in- 
volved. From previous measurements of the Whipple 
10 m electronics and the design of the SGARFACE 
DAQ, the expected digital-count (dc) per photo-electron 
(pe) is 0.92 ± 0.1, where the error is due to system- 
atic uncertainties of some 10 m measurements. As these 
measurements were done only once and at separate 
times, a second method was developed to track instru- 
mental variations. It consists of two parts: 

1) Events that are recorded simultaneously by both, 
Whipple and SGARFACE, are used to cross-calibrate 
the light-level between each cluster of pixels. Cosmic 
rays produce coincident events about 70% of the time, 
providing a large enough statistical sample to perform a 
nightly calibration. One such coincident event is shown 
in Fig. [23] The ratio of the total light measured in digital 
counts by Whipple versus SGARFACE is 3.54 ± 0.17, 
and it is very stable over the five-year time span. The 
statistical error is due to gain variations between pixels. 

2) The light-throughput of the Whipple 10 m is mea- 
sured from recorded images of complete muon rings. 
Cherenkov light detected from complete muon rings 
measures the quality of the local, < 500 m, atmo- 
sphere and the telescope optics and electronics through- 
put. This includes gain changes in the PMT The peri- 
odic variations in the measured light throughput seen in 
Fig.|24]is due to aging of the PMTs followed by subse- 



quent increase of the applied high voltage. Because of 
uncertainties in the ultra-violet response of the mirrors 
and PMTs, the dc/pe value was tied to the directly mea- 
sured value of 3.3±0.3 dc/pe for the 10 m in November 
2000 instead of using the absolute value derived from 
the muon calibration. 

Using this cross-calibration method, the dc/pe ratio 
was determined to be 0.93±0.1 dc/pe for SGARFACE 
at the time of the directed measurement and is therefore 
in agreement with the value obtained from the direct 
measurements. We note that two channels (18 and 54) 
exhibit large variations in the calibration ratio and we 
found this to be due to either, one, or two channels 
dropping out in the Whipple DAQ. 




Fig. 23. A cosmic-ray event detected simultaneously with SGAR- 
FACE (left) and Whipple (right). The total charge is represented by 
the size of the shaded area in each pixel. The extend of the fitted 
two-dimensional Gaussian is shown in each image by an ellipse. 
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Fig. 24. Light throughput measured relative to the calibrated point 
of 0.93 dc/pe in November 2000. 



6.2. Stability of the FADC Pedestals 

The SGARFACE FADC input is biased by pull-up 
resistors to avoid zero-crossing of the digital values. 
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This pedestal level is measured for every channel on 
an event-by-event basis from the first 750 samples with 
value below 40 dc. Typical pedestal values are 10-30 dc. 
Pedestal variations arise from the night-sky background. 
The pedestal level may also drift with temperature and 
noise in the data-acquisition system. 

The RMS pedestal fluctuation, as measured by the 
standard deviation of the distribution during an event, 
lies between 1.5 dc and04.5 dc, depending on night-sky 
condition and PMT gain. Changes in the mean pedestal 
value between consecutive events show a spread of 1.3 
dc RMS. Over the long term, the pedestals were very 
stable; all channels showed less than 2 dc drift over the 
five year analysis period. A minor hardware problem 
was identified where on rare occasions, traces would 
start at dc and remain at that level for up to entire 
stack length. The analysis filters out rare instances of 
missing FADC traces. 

6.3. Camera Pixel Gain 

The sensitivity across the field of view of the Whip- 
ple 10 m camera is set by the gain of individual PMTs. 
To achieve a uniform gain, the high-voltages of the 
PMTs are adjusted at the beginning of each observing 
season. In addition, the Whipple high-energy 7-ray ob- 
servations are flat-fielded in software with an auxiliary 
calibration run taken each night. For this purpose, a Ni- 
trogen flasher is used to uniformly illuminates the cam- 
era. These calibration data are not directly usable for 
analysis of SGARFACE data, because the flashes are 
too bright and saturate the FADCs. 

The gain uniformity of SGARFACE channels was 
tested with two methods: 

1) Dedicated flat-fielding data were taken at variable 
brightness levels using a laser system. The laser flashes 
were recorded by both, SGARFACE and the Whipple 
10 m DAQ. The flat-fielding coefficients were deter- 
mined independently for each system and are shown 
in Fig. [25] where the gains of the Whipple 10 m have 
been mapped onto the appropriate SGARFACE chan- 
nels. Both sets of flat-fielding coefficients are in agree- 
ment with each other. The x 2 -difference between the 
two distribution, calculated using the sample standard 
deviation of Whipple gain and the standard deviation of 
the mean SGARFACE gain, is 15.7 with 55 degrees of 
freedom. Therefore, the flat-fielding coefficients from a 
Whipple Nitrogen run can be used to predict the true 
flat-fielding coefficients of the SGARFACE pixel map. 
Figure |271 

The variation of the gain seen in Figf25]is shown in 



Fig. [26] in camera coordinates. As can be seen, during 
these flat-fielding runs the laser was not well aligned 
with the optical axis of the telescope. 

2) The Whipple 10 m flat-fielding coefficients as de- 
termined from the Nitrogen-flasher data are shown in 
Fig. [27] It has been established during the many years 
of operation of the 10 m, that the Nitrogen flasher uni- 
formly illuminates the field of view, so that any mea- 
sured gain differences are intrinsic to the camera. The 
distribution is uniform, except for two outliers, cor- 
responding to SGARFACE pixels 34 and 49. In both 
cases, the Whipple analysis reveals a low-gain PMT: 
PMT number 303 corresponding to SGARFACE pixel 
49 and PMT number 316 for pixel 34. 

Burst of 7 rays produce broad images of Cherenkov 
light that is spread over a large number of SGARFACE 
pixels. Typically, 10 to 30 SGARFACE pixels would 
contain light, with the gain for SGARFACE pixels vary- 
ing randomly between 0.8 and 1.2, see Fig. [27] With 
gains determined to within 20%, the parameter calcula- 
tion is effected very little since a large number of pixels 
are involved. For example, the total amount of detected 
light for a cluster of 15 pixels has an uncertainty of 
0.2/\/l5 = 0.05. 
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Fig. 25. Comparison of the flat-fielding coefficients derived with 
from SGARFACE and Whipple data. In the case of Whipple data, 
the individual pixels were remapped to the cluster making up the 
SGARFACE pixels. The data is from a one minute laser run on 
MJD 52821, Whipple run gt025056. Error bars show the standard 
deviation of the distribution. 
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